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Abstract 

We present the results of magnetic properties and entanglement for kagome 
lattice using Heisenberg model with two-, and three-site exchange interac- 
tions in strong magnetic field. Kagome lattice correspond to the third layer 
of fluid 3 He absorbed on the surface of graphite. The magnetic properties 
and concurrence as a measure of pairwise thermal entanglement are studied 
by means of variational mean-field like treatment based on Gibbs-Bogoliubov 
inequality. The system exhibits different magnetic behaviors, depending on 
the values of the exchange parameters (J2, J3). We have obtained the magne- 
tization plateaus at low temperatures. The central theme of the paper is the 
comparing the entanglement and magnetic behavior for kagome lattice. We 
have found that in the antiferromagnetic region behaviour of the concurrence 
coincides with the magnetization one. 

1 Introduction 

One can model solid and fluid 3 He films as the systems of almost localized identical 
fermions. Since the light mass spin-1/2 3 He atoms are subject to a weak attractive 
potential, the theoretical explanation of magnetism is based on the multiple-spin 
exchange mechanism pp. An important case is represented by solid and fluid 3 He 
films absorbed on the surface of graphite [21 El S] since it is a typical example of 
a two-dimensional frustrated quantum-spin system [5]. The first and second layers 
of the system form a triangular lattice [6], while the third one forms a system of 
quantum 1/2 spins on a kagome lattice [?]. 

The phenomenon of magnetization plateau has been studied during the past 
decade both experimentally and theoretically. The plateaus may be exhibited in the 
magnetization curves of quantum spin systems at very low temperatures in case of 
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strong external field. Magnetization plateaus appear in a wide range of models on 
chains, ladders, hierarchical lattices, theoretically analysed by dynamical, transfer 
matrix approaches and exact diagonalization in clusters (see Ref. [8]- [18]). In [T9] 
dynamical system theory has been used to study magnetization plateaus on the 
kagome chain with two-, three- and six-site exchange interactions. 

Recent years much effort has been put into studying the entanglement of mul- 
tipartite systems both qualitatively and quantitatively [201 [21] ■ Entanglement has 
gained renewed interest with the development of quantum information science. En- 
tangled states constitute a valuable resource in quantum information processing [22J. 
Numerous different methods of entanglement measuring have been proposed for its 
quantification [23]. In this paper we use concurrence [21] as entanglement measure 
of the spin-1/2 system. 

In the present paper mean-field like approach, based on the Gibbs-Bogoliubov 
inequality, was used to study entanglement and magnetic properties of a kagome 
lattice [25]. This method can also be applied to study thermal entanglement in 
many-body systems [26| 127] . 

The key result of the paper is concentrated on the comparison of specific (peaks 
and plateaus) features in magnetization and thermal entanglement properties in the 
above mentioned model using variational mean-field like Gibbs-Bogoliubov inequal- 



This paper is organized as follows: in Section 2 we introduce the Heisenberg 
model on the kagome lattice with two- and three-site exchange interactions. In 
Section 3 mean-field like approximation, based on Gibbs-Bogoliubov inequality, has 
applied on kagome lattice. The magnetic properties of the model are investigated 
in Section 4. In Section 5 concurrence as a measure of entanglement is studied and 
compared with magnetic properties of kagome lattice. The conclusive remarks are 
given in Sec. 6. 

2 Heisenberg model Hamiltonian on kagome lat- 
tice with two-, and three-site exchanges inter- 
actions 

The Hamiltonian of the Heisenberg model is 



where H ex represent spin exchanges and Hz is responsible for magnetism. The 
expression for Hz can be written as 



ity 



H = H ex + Hz, 



(1) 



H z 



(2) 
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Figure 1: Kagome lattice. 



where 7 ,B are gyromagnetic ratio and magnetic field. According to [T] the multiple 
spin exchanges Hamiltonian can be written as 

H ex = J 2 £ P M - J 3 E fa* + P ^lk) +■■■, (3) 

(i,j,k) 

where Pij,Pij,k represent the two-, and three-spin cyclic permutation operators. 
The sums are taken over all distinct two- and three-cycles. The expression of pair 
transposition operator P^ has been given by Dirac 

P,j = \( l + > ( 4 ) 

where er^ are the Pauli matrixes, acting on the spin at the z-th site. Using this 
expression for P^ one can find the expressions for the other operators of cyclic 
rearrangement, (see [UHS]) Here is the expression for three spin exchange operator: 



P i,j,k + P i,lk = „ (! + + (T .y CT A' + VkVi) ■ (5) 



As mentioned above the third layer of 3 He system is kagome lattice (see Fig 1). 
In kagome lattice each edge belongs to only one triangle and each site belongs to two 
triangles, therefore one can combine first two summations in ([3]) and Hamiltonian 
for kagome lattice can be written in the following form: 



Triangles . 



J2 — J3 / \ h 

(TitTj + <7,jCT k + (T k a-i) - -[a i + a i + a k 



(6) 



According to [6] the effective value of the exchange parameters J = J 2 — 2J3 on 
triangular lattice, which has been estimated experimentally from susceptibility and 
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specific-heat data for solid 3 He is J = —3.07mK. Therefore, the three-site exchange 
(J3) is dominant on the triangular lattice. This is a consequence of the fact that on 
triangular lattice, for high densities the probability of a triple permutations of 3 He 
atoms is dominant than a pair one. In the case of fluid 3 He, which is described by 
kagome lattice, pair exchanges become more probable, since every edge belongs to 
one triangle and one hexagon, in contrast to the triangular lattice for which every 
edge belongs to two triangles. Moreover, we did not take into account the six-site 
exchange interaction which is the antiferromagnetic one. Taking into account above 
mentioned facts we can consider effective pair exchange permutations (J 2) more 
dominant than three-site one (J3). 



3 Basic Gibbs-Bogoliubov mean-field formalism 

Here we apply the variational mean-field like treatment based on Gibbs-Bogoliubov 
inequality [25] to solve the Hamiltonian fl6]). This implies that the free energy 
(Helmholtz potential) of system is 

F<F + (H-H ) , (7) 

where H is the real Hamiltonian which describes the system and H is the trial 
one. F and Fq are free energies corresponding to H and Hq respectively and (. . .) 
denotes the thermal average over the ensemble defined by Hq. By introducing trial 
Hamiltonian for our model (Eq. ([6]) kagome lattice) containing unknown variational 
parameters one can minimize right hand side of Bogoliubov inequality (J7J) and get 
the values of those parameters, 

For antiferromagnetic interactions, the trial Hamiltonian will consist of two parts 
describing the two sublattices. We introduce a trial Hamiltonian Ho as a set of non- 
interacting clusters (triangles) on two sublattices in different external self-consistent 
fields: 



Ho = T, H o\ (8) 



where 



HP = A x (<j\<j\ + 0*0* + - lv x [(alY + (ai) z + (alY] , (9) 



where A and j v variational parameters, and Aj labels different noninteracting rect- 
angles (see Fig. 1, grey triangles) and 



7u = 7a for sublattice (a), 
7„ = 7 b for sublattice (b). 



(10) 



It should be emphasized that in trial Hamiltonian spins <j\, of the Aj-th triangle do 
not interact with the spins <7 3 k of the Aj triangle if i ^ j, therefore these spins are 
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statistically independent. Suppose the real Hamiltonian H (JS]) can be represented 
in the following form 

h = J2h^\ (ii) 

where ifW is the contribution of spins on the single triangle in real Hamiltonian 
H and index of summation Aj runs over the different triangles (see Fig. 1, grey 
triangles). Terms of real Hamiltonian <j\<j\ + cr\a\ + cr l 3 a\ must be included in 
but terms like cr l a cr^ (see Fig. 1 solid lines) should be included both in /fW and 
Consequently, has the following form: 



H® = h=h Q » + S ^M +S ^ +S ^ -h ±(KY, (12) 

Z \ r=2,3 Z fe=l,3 Z a=l,2 Z / r=l 

where 

a« = <r;a* + <r<<r* + <r|<ri. (13) 

In the expressions for we take half of each term er l a &1 because this term should 
be included in two different triangles. 

Inequality (j7]) can be rewritten now for the single triangle on each sublattice (v): 

fv <(fo) v + {H^-H^) o , (14) 

where H^> is the real and Hq the trial Hamiltonians of the triangle, f v and (fo) v free 
energies of the one triangle on sublattice (v) defined by and respectively. 
By denoting magnetizations of the sublattices (a) and (6) respectively m a and m h 
and taking into account that spins a l T belong to sublattice (a) and spins a^ k belong to 
sublattice (6) and fact that spins <r l T and erp (i ^ j, k) are statistically independent 
we get: ((4T' y ) = 0, m a = ((a*) 2 ) /2, m b = {{a^f) /2 and (<^) = ((<)*> x 
((^b) / = 2m a 2m;,. One can rewrite inequality (THl) as follows: 

/« < (fo)v + - Js ~ (a) + J2 - J3 6(2m a 2m fe ) - (h - 7^) 6m v . 
7 a = h - (J 2 - h)m b , 7b = h - (J 2 - h)m a . 
Minimizing the right hand side of (fT5|) in order to 7 a ,7& and A and using the fact, 

— = (a} and — 
parameters: 

A 



that -7^- = (a) and — — = — 6m„ we obtain the following values for the variational 



J2 — Ja 



2 ' 

7 a = h - (J 2 - J 3 )«2 b . 

76 = h - (J 2 - h)m a . (15) 

The Hamiltonian H$ was chosen to be exactly solved. By diagonalization Hamil- 
tonian Hq^ one can find eigenvectors end eigenvalues of the trial Hamiltonian 
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a) 



b) 



Figure 2: Magnetization m a versus external magnetic field h for J2 = 3 mK, J3 = 
2.5 mK. at a) T=0.15 mK b) T=0.01 mK. 

4 Magnetic properties 

Here and further exchange parameters (J2, J3) and magnetic field h is taking in 
Boltzman's constant scaling i.e. Boltzmann's constant is set to be ks = 1. 

The results of the previous section can be used for investigation of the magnetic 
properties of our model. The magnetization of arbitrary site is defined as 



m. 



Tr(S v e- H ' T ) 



(16) 



S v corresponding spin operator on sublattice (v), H is the Hamiltonian ([9]) with 
constants ( fl5l) and Z is partition function of the system. But according to ( fl5l) the 
Hamiltonian of sublattice (a) depends on m b through 7 a and vice versa. For defined 
above magnetization we obtain the following expression: 



m n 



1 3 sinh (^) + sinh (>-) + 2e(^)sinh fa* 
6 cosh (^) + cosh (^) + 2e(^)cosh (^) 



and 7 a = h-(J 2 -J3)"i6- (17) 



The dependance of magnetization m a from external magnetic field h can be found 
by solving the this recursive equation for each value of magnetic field h. 

At relatively high temperatures the recursive equation has one stable solution 
and therefore magnetization curves of sublattices (a) and (6) coincide (see Fig. 2(a)). 
With decreasing temperature the solution of recursive equation ceases to be stable 
and, therefore, the magnetization of different sublattices are no longer equal. The 
partially saturated phase emerges in form of the magnetization plateaus (see Fig. 
2(b) for T = 0.01 mK, J2 = 3 mK, J3 = 2.5 mK), which can be associated with 
a staggered magnetization or short range antiferromagnetism (AF) in frustrated 
kagome geometry. Indeed, the appearance of plateaus in magnetization curve at m = 
±1/6 can be explained as stability of trimeric states in available (ttl; tit; Itt) an d 
(t>Uo 4t4; til) configurations. Moreover, zero field magnetisation of one sublattice 
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m(h=0) 




Figure 3: Dependence of the magnetization m a (solid line) and concurrence (dashed 
line) at zero external field from temperature T at J 2 = 3 mK, J 3 = 2.5 mK. 

becomes nonzero. In figure 3 the solid line shows the temperature dependence of the 
magnetization in the absence of external magnetic field. The magnetization tents 
gradually to zero near the second-order transition temperature Tq between ordered 
m a 7^ and disordered m a = phases. 

5 Concurrence and thermal entanglement 

The mean-field like treatment transforms kagome lattice to the set of noninteracting 
triangles in effective field, therefore quantum correlations can be exactly accounted. 
This allows in terms of three-qubit XXX Heisenberg model in effective magnetic field 
7 to study the thermal entanglement properties. We will study the concurrence as a 
measure of pairwise entanglement [24J. The concurrence C(p) corresponding to the 
density matrix p is defined as 

C{p) = max{Ai - A 2 - A 3 - A 4 ,0}, (18) 

where Aj are the square roots of the eigenvalues of the operator 

~p = p l2 (ol®al)pl 2 {aX®al), (19) 

where p\ 2 = Tr%p is the reduced density matrix of the pair and p is defined in the 
following way 

z i=i 

where Z is the partition function of the system and ipi and Ei are eigenvectors and 
eigenvalues of the Hamiltonian Hq respectively (see Eq. [9]). p\ 2 has the following 
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a) b) 

Figure 4: Dependence for (a) magnetization m and (b) concurrence C(p) versus the 
magnetic field h and the coupling constant J 2 at J3 = 2.5 mK and T = 0.2 mK. 



form 

P12 



( u ^ 

w y 

y iu 

V v ) 



(21) 



where u,w,y and u are some functions of variables 7, A and T. Using ( fl8|) . ( fl9|) and 
( |2~TT) one can find the following expression for the concurrence C(p) 

C{p) = max{\y\ - V^v, 0}. (22) 

In this equation one must replace 7 with h — 2(J 2 — Js)m (see Eq.( fT5|) ). therefore the 
concurrence C(p) is the function of magnetisation m. To calculate concurrence one 
must solve transcendental equation (fT7|) for each set of parameter values ( J 2 , J3, h, T) 
and insert corresponding solution to the equation (!22l . 

It is curious to discuss some similarities of statistical and quantum characteristics 
of our system. We consider magnetization as a statistical characteristic. In figure 
3(a) plotted the magnetization as a function of the coupling constant J 2 (for fixed 
value of J3 = 2.5 mK) and the external field h, at a relatively high temperature 
T = 0.2 mK. As a quantum characteristic we consider entanglement (concurrence 
C(p)). In figure 3(b) the concurrence as a function of the J 2 (J3 = 2.5 mK) is 
shown for the same value of temperature. Our calculations show that the magnetic 
characteristics is similar to that of bipartite entanglement. Indeed, comparison 
of figures 3(a) and 3(b) shows that regions corresponding to the magnetization 
plateaus, coincide with the plateaus on concurrence plot. 
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6 Conclusion 



In this paper we find strong correlations between magnetic properties and quantum 
entanglement in the Heisenberg model with two-, and three-site exchange interac- 
tions in strong magnetic field on the kagome lattice, which correspond to the third 
layer of fluid 3 He absorbed on the surface of graphite. We adopted variational 
mean-field-like treatment (based on the Gibbs-Bogoliubov inequality) of separate 
clusters in effective magnetic fields and studied magnetic properties and concur- 
rence as a measure of pairwise thermal entanglement. The system exhibits different 
magnetic behaviors, depending on the values of the exchange parameters {J2,J%)- 
We have obtained the magnetization plateaus at low temperatures. We have found, 
that in the antiferromagnetic region behavior of the concurrence coincides with the 
magnetization one. The comparison of magnetization and concurrence shows that 
regions corresponding to the magnetization plateaus, coincide with the plateaus on 
concurrence plot. 
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